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1.  Introduction 


Conventional  transparent  armor  consists  of  glass  laminates  with  polymer  interlayer  and  backing. 
It  has  been  demonstrated  that  the  materials,  the  ratio  of  materials,  and  the  type  and  thickness  of 
the  interlayers  affect  the  ballistic  efficiency  of  the  laminate. 

Borosilicate  glass,  Starphire*  ultra-clear  soda-lime  glass,  and  the  transparent,  polycrystalline 
ceramic  AlON  are  materials  being  considered  for  transparent  armor  applications.  A 
comprehensive  series  of  Edge-on  Impact  (EOI)  tests  has  been  conducted  in  order  to  examine 
wave  and  damage  propagation  through  the  single  materials  (baseline  tests)  and  laminated 
structures. 

The  numerical  simulation  part  of  this  project  focused  on  the  modeling  and  simulation  of 
projectile  impact  on  the  poly  crystalline,  transparent  ceramic  AlON.  In  the  approach  chosen  here, 
the  ceramic  was  not  treated  as  a  continuum.  In  order  to  model  the  polyhedral  microstructure  of 
the  ceramics,  a  code,  which  is  based  on  the  theory  of  power  diagrams,  was  implemented.  After 
having  obtained  a  three-dimensional  (3-D)  grain  structure  that  corresponds  on  average  to  what  is 
observed  in  micrographs,  these  structures  are  being  meshed  in  3-Ds  using  tetrahedra  in  3-D  vs. 
triangles  in  two  dimensions.  The  generated  and  meshed  microstructures  are  used  as  input  for  a 
commercial  program  package  (LSDyna)  and  subsequent  finite-element  (FE)  analyses. 


2.  Experimental  Set-up 


The  EOI  test  method  coupled  with  a  high-speed  Cranz-Schardin  camera,  with  frame  rates  up  to 
107  fps,  has  been  applied  in  order  to  visualize  damage  propagation  and  dynamic  fracture  in 
structural  ceramics.  Two  different  optical  configurations  were  employed.  A  regular  transmitted 
light  shadowgraph  set-up  was  used  to  observe  wave  and  damage  propagation  and  a  modified 
configuration,  where  the  specimens  were  placed  between  crossed  polarizers  and  the  photo-elastic 
effect  was  utilized  to  visualize  the  stress  waves.  Pairs  of  impact  tests  at  approximately 
equivalent  velocities  were  carried  out  in  transmitted  plane  (shadowgraphs)  and  crossed  polarized 
light.  Figure  1  shows  a  schematic  of  the  EOI  test  with  the  added  crossed  polarizers;  figure  2 
illustrates  an  exploded  view  of  the  impactor/sample  interaction. 


* 

Starphire  is  a  registered  trademark  of  PPG,  Pittsburgh,  PA. 
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Figure  1.  EOI  test  set-up  with  Cranz-Schardin  camera. 


Figure  2.  Close-up  view  of  test  sample  set-up  for 
shadowgraphs. 


Both  steel  solid  cylinder  and  spherical  impactors  have  been  used  at  velocities  from  350  to  450  m/s 
on  100-  x  100-  x  10-mm  plates.  Once  the  baseline  glass  materials  were  tested  and  analyzed,  multi- 
component  glass  laminates  were  produced  and  tested  at  ~400  m/s.  The  data  collected  from  the 
EOI  test  consists  of  a  series  of  20  photographs  as  a  function  of  time,  typically  at  0.25-2  ps 
intervals.  Pairs  of  impact  tests  at  approximately  equivalent  velocities  are  carried  out  in  plane  and 
crossed  polarized  light  to  correlate  the  dynamic  fracture  with  the  associated  stress  fields. 

Detailed  graphs  are  then  created  plotting  crack,  damage,  and  compression  and  shear  stress  wave 
velocities.  More  details  on  the  set-up  can  be  found  in  the  first  interim  report  (7). 
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3.  Nomenclature 


Since  two  different  methods  (shadowgraphs,  crossed  polarizers)  were  employed  for  the 
visualization  of  wave  and  damage  propagation,  where  not  only  the  stress  waves  appear  in 
different  forms,  but  also  different  types  of  fracture  can  be  distinguished,  it  is  necessary  to  clarify 
the  designations  of  the  various  phenomena. 

3.1  Fracture 

The  nomenclature  for  the  different  fracture  types  and  patterns  was  defined  in  earlier  EOI  studies 
(2)  and  is  also  used  in  the  work  reported  here.  Figures  3  and  4  illustrate  typical  damage  patterns 
observed  with  soda-lime  glass  at  low  (—50  m/s)  and  medium  (~200  m/s)  impact  velocities. 


Secondary  cracks 


Primary  cracks 


Crack  center 
due  to  nucleation 


Figure  3.  Typical  damage  pattern  in  float  glass  at  vP  ~  50  m/s; 

cR  =  rayleigh  wave  velocity;  cT  =  transversal  wave  velocity; 
vCr  =  terminal  crack  velocity. 
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Figure  4.  Typical  damage  pattern  in  float  glass  at  vP  ~  200  m/s. 


While  figure  5  illustrates  the  typical  damage  pattern  in  Starphire  high-purity  soda-lime  glass  at 
an  impact  velocity  of  ~400  m/s,  figure  6  shows  the  damage  pattern  in  AlON,  transparent, 
polycrystalline  ceramic  at  about  the  same  impact  velocity. 


Figure  5.  Typical  damage  pattern  with  Starphire  glass  at 
vP  ~  400  m/s. 
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Figure  6.  Typical  damage  pattern  in  A10N  at  vP  ~  400  m/s. 


3.2  Waves 

The  two  different  optical  techniques  employed  exhibit  different  sensitivities  with  respect  to  the 
stress  level  that  can  be  visualized.  In  a  shadowgraph  image,  the  light  intensity  depends  on  the 
second  spatial  derivative  d2n /dx2  of  the  refractive  index,  whereas  in  the  crossed  polarizers  set-up, 
the  intensity  of  the  transmitted  light  depends  on  the  photo-elastic  properties  of  the  material. 
Therefore,  it  is  possible  that  the  first  visible  wave  front  in  the  shadowgraph  configuration 
appears  at  a  different  position  than  the  forefront  of  the  stress  wave,  visible  in  the  crossed 
polarizers  set-up.  Both  techniques  can  visualize  different  parts  of  the  same  stress  wave.  This  is 
demonstrated  in  figure  7,  which  shows  the  shadowgraph  and  the  corresponding  crossed 
polarizers  photograph  of  a  Starphire  specimen  at  8.7  ps  after  impact. 

The  width  of  the  stress  wave  can  be  estimated  from  the  length  of  the  projectile  and  the 
longitudinal  wave  speed  in  the  projectile  and  target  material.  When  the  projectile  hits  the  target, 
a  wave  is  generated  not  only  in  the  target  but  also  in  the  projectile.  The  maximum  length  of  the 
stress  pulse  is  given  by  the  time  the  wave  needs  to  travel  once  through  the  projectile  and  back. 
The  thickness  of  the  solid  cylindrical  part  of  the  projectile  was  8  mm.  With  a  longitudinal  wave 
speed  of  5100  m/s  in  steel,  this  yields  a  stress  pulse  length  of 
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Figure  7.  Shadowgraph  (top)  and  corresponding  crossed 
polarizers  photograph  (bottom)  of  Starphire 
specimen,  8.7  jus  after  impact  at  ~  400  m/s. 


Release  waves  from  the  edges  of  the  projectile  are  neglected  in  this  estimate.  The  longitudinal 
wave  velocity  in  Starphire  glass  is  -5800  m/s.  Therefore,  the  maximum  spatial  width  of  the 
stress  wave  in  the  target  can  be 


As 


stress 


„  ,  _ mm  , 

~3.1  ps-5.8 - ~18  mm. 

ps 


(2) 
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The  formation  of  the  second  wave  front  is  caused  by  the  geometry  of  the  specimens  and  is  due  to 
the  initiation  and  superposition  of  waves,  generated  at  the  side  surfaces  of  the  specimens.  The 
distance  to  the  first  wave  front  depends  on  the  thickness  and  Poisson’s  ratio  (2). 


4.  Baseline  Results 


Experiments  performed  in  plane  light  show  the  evolution  of  damage  and  material  failure,  while 
the  photoelastic  visualization  illustrates  the  stress  wave  propagation  as  a  function  of  time. 

Figure  8  shows  a  selection  of  two  shadowgraphs  (top)  and  corresponding  crossed  polarizers 
photographs  (bottom)  of  a  baseline  test  with  Starphire  glass,  impacted  by  a  spherical  steel 
projectile  with  a  16  mm  diameter  at  440  m/s.  The  shadowgraphs  show  a  crack  front  growing 
from  the  impacted  edge  of  the  specimen.  Only  one  crack  center  can  be  observed  close  to  the 
upper  edge  of  the  specimen.  The  crossed  polarizers  photographs  illustrate  the  propagation  of  the 
longitudinal  and  the  transversal  stress  waves.  Release  waves  due  to  reflections  at  the  upper  and 
lower  edge  can  also  be  recognized.  Note  that  damage  appears  dark  on  the  shadowgraphs  and  the 
zones  with  stress  birefringence  are  exhibited  as  bright  zones  in  the  crossed  polarizers 
photographs. 


Figure  8.  Selection  of  two  shadowgraphs  (top)  and  crossed  polarizers  photographs 
(bottom)  from  impact  on  Starphire  glass  with  steel  sphere  at  440  m/s. 
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Figure  9  shows  a  selection  of  two  shadowgraphs  along  with  the  corresponding  crossed  polarizers 
photographs  of  the  baseline  tests  with  the  cylindrical  projectile.  A  coherent  damage  zone  is 
growing  from  the  impacted  edge,  preceded  by  a  zone  with  separated  crack  centers,  initiated  by 
the  stress  waves.  It  can  be  recognized  that  the  stress  wave  front  appears  more  advanced  and 
exhibits  a  different  curvature  in  the  crossed  polarizers  view.  This  seeming  discrepancy  can  be 
explained  by  the  different  sensitivities  that  the  different  optical  techniques  employed  exhibit  with 
respect  to  the  stress  level  that  can  be  visualized.  In  a  shadowgraph  image,  the  light  intensity 
depends  on  the  second  spatial  derivative  d2n/dx2  of  the  refractive  index,  whereas  in  the  crossed 
polarizers  set-up,  the  intensity  of  the  transmitted  light  depends  on  the  photo-elastic  properties  of 
the  material.  Therefore,  it  is  possible  that  the  first  visible  wave  front  in  the  shadowgraph 
configuration  appears  at  a  different  position  than  the  forefront  of  the  stress  wave,  visible  in  the 
crossed  polarizers  set-up.  Both  techniques  can  visualize  different  parts  of  the  same  stress  wave. 


Figure  9.  Selection  of  two  shadowgraphs  (top)  and  crossed  polarizers 

photographs  (bottom)  from  impact  on  Starphire  glass  with  steel 
cylinder  at  390  m/s. 

The  high-speed  photographs  from  baseline  tests  with  the  steel  sphere  and  cylinder  on  Boro  float* 
glass  are  presented  in  figures  1 0  and  1 1 .  Crack  and  damage  patterns  very  similar  to  those 
observed  with  soda-lime  glass  can  be  recognized.  However,  the  crack  and  damage  velocities 
determined  from  the  position-time  analysis  of  the  high-speed  photographs  were  higher  than  with 


Borofloat  is  a  registered  trademark  of  Schott  North  America,  Inc.,  Elmsford,  NY. 
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Starphire  glass.  Tables  1  and  2  summarize  the  crack  velocity,  damage  velocity,  and  wave 
propagation  data  of  the  baseline  tests  with  Starphire  and  Borofloat  glass,  respectively. 


Figure  10.  Selection  of  four  shadowgraphs  from  impact  on 
Borofloat  glass  with  steel  sphere  at  430  m/s. 


Figure  11.  Selection  of  two  shadowgraphs  (top)  and  crossed 
polarizers  photographs  (bottom)  from  impact  on 
Borofloat  glass  with  steel  cylinder  at  390  m/s. 
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Table  1.  Measured  wave,  crack,  and  damage  velocities  with  Starphire  soda-lime  glass. 


Impactor 

Optical 

Set-up 

Long.  Wave 
Velocity  cL 
(m/s) 

Trans.  Wave 
Velocity  cT 
(m/s) 

Crack 
Velocity  vCr 
(m/s) 

Damage 
Velocity  vD 
(m/s) 

Sphere 

Shadowgraph 

— 

— 

1580 

1580 

Sphere 

Crossed  polarizers 

5763 

3518 

— 

— 

Cylinder 

Shadowgraph 

5761 

— 

— 

3270 

Cylinder 

Crossed  polarizers 

5779 

— 

— 

— 

Table  2.  Measured  wave,  crack,  and  damage  velocities  with  Borofloat  glass. 


Impactor 

Optical 

Set-up 

Long.  Wave 
Velocity  cL 
(m/s) 

Trans.  Wave 
Velocity  cT 
(m/s) 

Crack 
Velocity  vCr 
(m/s) 

Damage 
Velocity  vD 
(m/s) 

Sphere 

Shadowgraph 

5462 

— 

2034 

2034 

Cylinder 

Shadowgraph 

5531 

— 

— 

4150 

Cylinder 

Crossed  polarizers 

5635 

— 

— 

— 

5.  Starphire  Glass  Laminates 


The  influence  of  a  polyurethane  (PU)  bonding  layer  on  wave  and  damage  propagation  was 
examined  with  cylindrical  projectiles  only.  Four  pairs  of  tests  with  specimens  consisting  of  two 
parts  of  the  dimensions  50  x  100  x  9.5  mm  were  conducted  in  order  to  examine  the  influence  of 
interlayer  thickness.  Starphire  specimens  with  bonding  layers  of  thickness  0.64,  1.27,  2.54,  and 
5.08  mm  were  examined.  The  influence  of  two  PU  bonding  layers  was  tested  with  specimens 
that  were  built  of  three  parts  of  the  dimensions  30  x  100  x  9.5  mm. 

One  pair  of  tests  was  conducted  with  a  different  type  of  PU  interlayer,  designated  IM  500. 

Table  3  shows  a  compilation  of  these  tests. 

The  phenomenology  of  damage  and  wave  propagation  in  Starphire  laminates  has  already  been 
discussed  in  the  first  and  second  Interim  Report  (7,  3).  Therefore,  a  summary  of  these  results 
will  be  given  in  the  following  section.  Only  the  results  with  the  IM  500  type  bonding  layer  will 
be  shown  in  detail. 

Figure  12  illustrates  a  comparison  of  wave  propagation  and  damage  in  Starphire  specimens  with 
bonding  layers  of  thickness  0.64,  2.54,  and  5.08  mm.  The  impact  velocity  was  380  ±  5  m/s  in  all 
tests.  The  upper  line  of  pictures  shows  the  shadowgraphs,  while  the  corresponding  crossed 
polarizers  photographs  are  presented  in  the  lower  line  of  pictures,  respectively.  Figure  12a 
illustrates  the  specimens  at  10.7  ps  and  figure  12b  at  23.7  ps  after  projectile  impact.  The 
shadowgraphs  at  10.7  ps  show  that  the  first  glass  layer  (left  part  of  specimen)  is  damaged 
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Table  3.  Tests  with  laminated  Starphire  glass. 


Test  No. 

Interlayer 

(in) 

Impact  Velocity 

(m/s) 

Optical  Set-up 

Time  Intervals 

Picture  No.:  At 
(fis) 

15274 

0.025  DF 

375 

Shadowgraph 

1-14;  1;  14-20:  2  jus 

15275 

0.025  DF 

387 

Crossed  polarizers 

1-14:  1;  14-20:  2  ps 

15299 

0.050  DF 

389 

Shadowgraph 

1-14;  1;  14-20:  2  jus 

15300 

0.050  DF 

395 

Crossed  polarizers 

1-14:  1;  14-20:  2  ps 

15273 

0.100  DF 

384 

Shadowgraph 

1-16:  1;  16-20:  2  jus 

15280 

0.100  DF 

382 

Crossed  polarizers 

1-16:  1;  16-20:  2  ps 

15278 

0.200  DF 

385 

Shadowgraph 

1-10:  1;  10-20:  2  jus 

15279 

0.200  DF 

384 

Crossed  polarizers 

1-10:  1;  10-20:  2  ps 

15309 

2x0.100  DF 

404 

Shadowgraph 

1-8:2;  8-16:  1;  16-20:2 

15310 

2x0.100  DF 

395 

Crossed  polarizers 

1-8:2;  8-16:  1;  16-20:2 

15307 

0.100  IM  500 

388 

Shadowgraph 

1-4:  2;  4-8:  1;  8-20:  2 

15308 

0.100  IM  500 

399 

Crossed  polarizers 

1-4:  2;  4-8:  1;  8-20:  2 

0.64  mm 
1 - 


2.54  mm  5.08  mm 


0.64  mm 


2.54  mm 


5.08  mm 


(a)  10.7  jlxs 


(b)  23.7  jlis 


Figure  12.  Starphire  laminates  with  interlayer  of  different  thickness  impacted  by  steel  cylinder  at  380  m/s. 


through  the  coherent  fracture  front  growing  from  the  impacted  edge  and  through  the  nucleation 
of  crack  centers,  initiated  by  the  longitudinal  stress  waves.  At  that  time,  no  damage  can  be 
recognized  in  the  second  glass  layer  (right  part  of  specimen).  The  crossed  polarizers 
photographs  demonstrate  that  the  first  longitudinal  stress  pulse  has  not  yet  crossed  the  thickest 
glue  interlayer  (right),  whereas  the  stress  wave  is  clearly  visible  in  the  right  half  of  the  specimens 
with  the  thinner  glue  interlayer. 

After  23.7  ps  (figure  7b),  the  compressive  stress  pulse  has  already  been  reflected  as  a  tensile 
wave  at  the  rear  edge  of  the  specimens  in  all  three  cases.  The  shadowgraphs  illustrate  that 
damage  in  the  second  glass  layer  is  mainly  due  to  the  tensile  wave  and  starts  from  the  rear  edge 
of  the  specimen.  In  the  case  of  the  thickest  glue  interlayer,  only  little  damage  was  observed  in 
the  second  glass  layer. 
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The  effect  of  two  bonding  layers  of  2.54  mm  thickness  is  demonstrated  in  figure  13  which  shows 
a  selection  of  four  shadowgraphs  and  corresponding  crossed  polarizers  photographs  in  the  time 
interval  from  6  to  25  ps  after  impact  of  a  steel  cylinder  at  about  400  m/s.  The  first  layer  of  glass 
was  completely  damaged  within  the  first  15  ps.  Damage  could  be  recognized  in  the  second  layer 
after  16  ps,  when  the  first  crack  centers  became  visible  which  were  initiated  by  the  reflection  of 
the  compression  wave  at  the  interface  between  the  second  glass  and  the  second  bonding  layer. 

No  damage  was  observed  in  the  third  glass  layer  during  the  time  interval  of  observation. 


Figure  13.  Selection  of  four  shadowgraphs  and  corresponding  crossed  polarizers  photographs; 
Starphire  specimen  with  two  PU  bonding  layers;  impact  velocity  404/395  m/s. 


The  wave  propagation  in  the  specimens  was  analyzed  and  the  path-time  history  for  the  case  with 
the  two  2.54-mm  bonding  layers  is  presented  in  figure  14. 

The  results  with  an  aliphatic  PU  bonding  layer  of  type  IM  500  and  a  nominal  thickness  of 
0.64  mm  are  shown  in  figure  15  (selection  of  high-speed  photographs)  and  figure  16  (path-time 
history  of  wave  propagation). 

The  high-speed  photographs  in  figure  15  illustrate  that  the  damage  formation  in  the  second  glass 
layer  is  similar  to  the  case  with  the  DF  polyurethane  bonding  layer  of  the  same  thickness. 

When  the  waves  hit  an  interlayer,  one  part  is  reflected  while  the  other  part  is  transmitted  into  the 
next  glass  layer.  Due  to  the  low  acoustic  impedance  of  the  interlayer  compared  to  the  glass,  the 
amplitude  of  the  stress  pulses  is  attenuated  considerably.  The  low  wave  velocity  in  the 
interlayers  effected  a  time  delay  at  each  bonding  layer  compared  to  the  unperturbed  propagation 
through  the  glass. 
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Figure  14.  Path-time  history  of  wave  propagation  in  specimens  with  two  2.54-mm 
polyurethane  (DF)  bonding  layers. 


Figure  15.  Selection  of  four  shadowgraphs  and  corresponding  crossed  polarizers  photographs; 
Starphire  specimen  IM  500  bonding  layer;  impact  velocity  388/399  m/s. 
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Figure  16.  Path-time  history  of  wave  propagation  in  specimens  with  0.64-mm  IM  500  bonding 
layer. 

The  delay  times  measured  in  all  tests  were  plotted  in  a  delay  time  versus  bonding  layer  thickness 
diagram  (figure  17).  Linear  regression  of  the  data  yielded  an  average  delay  time  of  0.33  ps/mm. 
This  is  in  good  agreement  with  the  calculated  value  based  on  a  longitudinal  wave  velocity 
Cl  =  5770  m/s  for  Starphire  glass  and  Cl  ~  2000  m/s  for  the  polyurethane  ( 4 ). 
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6.  Influence  of  Interface  Shape 


A  series  of  EOI-tests  was  performed  in  order  to  study  the  influence  of  the  shape  of  the  interface 
between  Starphire  glass  layers  on  wave  and  damage  propagation.  Three  different  types  of 
interface  were  considered:  saw  tooth,  corrugated,  and  wave  shape.  The  different  interface  types 
are  illustrated  in  figure  18.  The  contour  of  the  wave  shape  interface  was  adopted  from  the 
contours  observed  in  the  baseline  tests  (crossed  polarizers)  with  Starphire  glass.  One  additional 
baseline  test  with  a  straight  interface  was  conducted  with  a  Starphire  specimen  of  12.3-mm 
thickness.  Table  4  shows  a  compilation  of  experimental  parameters  from  this  test  series. 


Figure  18.  Interface  types  tested. 


Table  4.  Compilation  of  experimental  parameters  from  tests  with  differently  shaped  interfaces. 


Test  No. 

Interface 

Imp. 

Velocity 

(m/s) 

Optical 

Set-up 

Time  intervals 

Picture  No.:  At 
(ps) 

Type 

DF  Thickness 
(mm) 

15724 

Saw  tooth 

3.9 

400 

Shadowgraph 

1-14:  1  jlis;  14-20:  2  jus 

15725 

Saw  tooth 

3.9 

342 

Crossed  polarizers 

1-14:  1  jlis;  14-20:  2  jlis 

15727 

Corrugated 

1.3 

405 

Shadowgraph 

1-14:  1  ps;  14-20:  2  ps 

15728 

Corrugated 

1.3 

403 

Shadowgraph 

1-14:  1  jlis;  14-20:  2  jlis 

15729 

Straight 

3.5 

404 

Shadowgraph 

1-14:  1  jlis;  14-20:  2  jlis 

15730 

Wave  shape 

1.8 

384 

Crossed  polarizers 

1-14:  1  jlis;  14-20:  2  jlis 

15731 

Wave  shape 

1.8 

404 

Shadowgraph 

1-7:  2  jlis;  7-17:  0.5  jlis 
17-20:  1  jlis 

Figure  19  shows  a  selection  of  four  shadowgraphs  along  with  the  corresponding  crossed 
polarizers  photographs  of  the  tests  with  the  saw  tooth  interface.  The  complete  sets  of  high-speed 
photographs  are  presented  in  figures  A-l  and  A-2  in  the  appendix. 
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Figure  19.  Selection  of  four  shadowgraphs  and  corresponding  crossed  polarizers  photographs; 
Starphire  specimens  with  saw  tooth  interface;  impact  velocities  400/342  m/s. 


In  contrast  to  the  results  of  all  tests  with  straight  interfaces  damage  by  nucleation  of  crack  centers 
can  be  observed  in  the  second  layer,  before  the  stress  wave  is  reflected  at  the  rear  edge.  Shortly 
after  the  stress  wave  has  crossed  the  interface,  cracks  can  be  recognized,  initiated  especially  at 
the  tips  of  the  saw  tooth  interface.  The  path-time  history  of  wave  and  damage  propagation  is 
depicted  in  figure  20.  The  best  fit  straight  line  (dashed)  through  the  nucleation  sites  of  three 
leading  crack  centers  yields  a  damage  velocity  of  4256  m/s  in  the  second  glass  layer. 


0  5  10  15  20  25 

t(ps) 


Figure  20.  Path-time  history  of  wave  propagation  in  specimens  with  saw  tooth  shaped 
interface. 
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Figure  21  shows  a  selection  of  four  shadowgraphs  along  with  the  corresponding  crossed 
polarizers  photographs  of  the  tests  with  the  corrugated  interface.  The  complete  sets  of  high¬ 
speed  photographs  are  presented  in  figures  A-3  and  A-4  in  the  appendix. 


Figure  21.  Selection  of  four  shadowgraphs  and  corresponding  crossed  polarizers  photographs; 
Starphire  specimens  with  corrugated  interface;  impact  velocities  405/403  m/s. 


The  sequence  of  high-speed  photographs  shows  that  cracks  are  initiated  all  along  the  interface, 
which  form  a  dense  crack  front,  propagating  at  the  terminal  crack  speed  of  1524  m/s  in  Starphire 
glass  (see  figure  22).  The  crack  front  in  the  second  glass  layer  is  preceded  by  several  crack 
centers,  nucleated  by  the  stress  wave,  before  it  is  reflected  at  the  rear  edge  of  the  specimen.  The 
path-time  history  of  wave  and  damage  propagation  is  depicted  in  figure  22. 


Figure  22.  Path-time  history  of  wave  propagation  in  specimens  with  corrugated  interface. 
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A  selection  of  four  shadowgraphs  along  with  the  corresponding  crossed  polarizers  photographs 
of  the  tests  with  the  wave  shaped  interface  is  depicted  in  figure  23.  Note  that  the  time  intervals 
chosen  in  the  test  with  the  crossed  polarizers  were  different  from  those  in  the  test  with  the 
shadowgraph  arrangement.  The  complete  sets  of  high-speed  photographs  are  presented  in 
figures  A-5  and  A-6  in  the  appendix. 


Figure  23.  Selection  of  four  shadowgraphs  and  corresponding  crossed  polarizers 

photographs;  Starphire  specimens  with  wave  shaped  interface;  impact  velocities 
384/404  m/s. 

As  observed  with  the  corrugated  interface,  cracks  were  initiated  all  along  the  wave  shaped 
interface,  which  formed  a  dense  crack  front,  propagating  at  a  speed  of  1514  m/s  (see  path-time 
history,  figure  24). 

The  sets  of  high-speed  photographs  from  the  tests  with  modified  interfaces  have  demonstrated 
that  damage  formation  occurred  in  the  second  layer  before  the  stress  waves  were  reflected  at  the 
rear  edge  of  the  specimens.  This  means  that  the  modified  interfaces  were  less  effective  with 
respect  to  damage  suppression  or  retardation  compared  to  the  laminates  with  straight  interfaces. 
In  contrast  to  the  straight  interfaces,  the  glass  surfaces  were  not  polished  with  the  modified 
interfaces.  Therefore,  initial  defects  at  the  surfaces,  due  to  machining  of  the  specimens,  could 
have  caused  stress  concentrations  and  crack  formation.  The  crack  initiation  at  the  tips  of  the  saw 
tooth  interface  can  possibly  be  attributed  to  this  mechanism.  However,  the  surfaces  of  the  one 
thick  specimen  with  straight  interface,  which  was  used  as  a  reference,  were  not  polished  either 
and  only  little  crack  formation  was  observed  in  the  second  layer  before  reflection  of  the  stress 
wave  at  the  rear  edge.  Figure  25  shows  a  selection  of  four  shadowgraphs  from  this  baseline  test. 
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Figure  24.  Path-time  history  of  wave  propagation  in  specimens  with  wave  shaped  interface. 


Figure  25.  Selection  of  four  shadowgraphs  from  baseline  test  with  thick  Starphire  specimen 
with  straight  interface;  impact  velocity  404  m/s. 


A  specimen  with  two  straight  DF  PU  interfaces  of  1 .24  mm  was  impacted  in  such  a  way  that  the 
interfaces  were  parallel  to  the  shot  axis  and  the  projectile  only  hit  the  intermediate  glass  layer 
(see  figure  26).  This  impact  configuration  mimics  the  situation  when  a  projectile  hits  one  tile  of 
a  mosaic  of  ceramic  tiles.  The  EOI  technique  allows  for  the  observation  of  the  wave  and  damage 
propagation  in  the  adjacent  tiles.  Figure  27  shows  the  complete  set  of  high-speed  photographs, 
which  were  taken  with  the  crossed  polarizers  set-up. 


Figure  26.  Schematic  of  EOI  configuration  to  mimic 
impact  on  one  tile  of  a  ceramic  mosaic. 
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Figure  27.  Complete  set  of  high-speed  photographs  in  crossed  polarizers  arrangement  from  EOI  on 
intermediate  glass  layer,  test  no.  15726. 
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7.  Results  With  AlON 


Two  pairs  of  tests  were  conducted  with  AlON  in  order  to  complement  the  earlier  test  series  (5). 
The  test  parameters  are  listed  in  table  5. 

Table  5.  Tests  with  AlON. 


Test  No. 

Projectile 

Test  Type 

Imp. 

Velocity 

(m/s) 

Opt. 

Set-up 

Time  intervals 

Picture  No.:  At 
(ps) 

15281 

Cylinder 

Confinement 

399 

Shadowgraph 

1-5:  1;  5-20:  0.5 

15283 

Cylinder 

Confinement 

394 

Crossed  polarizers 

1-5:  1;  5-20:  0.5 

15295 

Sphere 

Baseline 

429 

Shadowgraph 

1-5:  1;  5-20:  0.5 

15296 

Sphere 

Baseline 

426 

Crossed  polarizers 

1-5:  1;  5-20:  0.5 

In  the  first  pair  of  tests,  the  influence  of  a  steel  confinement  was  examined.  Since  the  acoustic 
impedances  of  steel  Tsteei  =  7.85T03kg/m3-5100  m/s  =  40T06  kg/m2s  and  AlON  Taion  = 
3.69T03kg/m3T0300  m/s  =  38-106  kg/m2s  nearly  match,  a  significant  reduction  of  the  amplitude 
of  the  reflected  waves  was  expected,  and  therefore,  a  reduced  influence  of  the  reflected  waves  on 
the  damage  pattern.  For  this  purpose,  a  target  mount  was  designed,  with  an  opening  only  for  the 
projectile.  Figure  28  shows  a  photograph  of  the  target  mount  with  an  AlON  specimen  alone,  and 
with  the  AlON  specimen  between  crossed  sheet  polarizers.  The  surfaces  of  the  square  bar  steels 
that  were  in  contact  with  the  AlON  specimen  were  ground  in  order  to  improve  the  contact 
between  the  two  materials.  The  photograph  with  specimen  between  crossed  polarizers  illustrates, 
that  already  relatively  low  stresses,  caused  by  the  mounting  of  the  target,  could  be  visualized. 


Figure  28.  Photographs  of  target  mount  with  AlON  specimen  (left)  and  mounted  specimen 
between  crossed  (sheet)  polarizers  (right). 
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In  the  first  test  with  confinement  the  high-speed  camera  was  triggered  too  late,  so  that  only  a  part 
of  the  sequence  of  photographs  showed  the  damage  development  in  the  specimen.  However, 
from  the  photographs  of  the  test  with  crossed  polarizers,  not  only  the  wave  front  could  be 
recognized  but  also  the  development  of  the  fracture  front  could  be  analyzed. 

A  selection  of  eight  photographs  from  the  test  with  steel  confinement  in  the  crossed  polarizers 
configuration  is  presented  in  figure  29. 


Figure  29.  Selection  of  eight  crossed  polarizers  photographs  from  test  no.  15283  with  steel 
confinement  and  impact  velocity  of  394  m/s. 


The  photographs  show  the  nearly  semicircular  shaped  stress  wave  (bright),  followed  by  a 
coherent  fracture  front  (dark),  which  is  preceded  by  single  crack  centers.  A  clear  boundary 
between  the  coherent  fracture  zone  and  the  damage  nucleation  zone,  where  the  crack  centers 
appear,  can  not  be  recognized.  The  only  significant  difference  in  the  total  damage  pattern 
compared  to  the  tests  without  confinement  is  the  growth  of  cracks  from  the  upper  and  lower  edge 
of  the  specimen.  This  crack  growth  could  be  due  to  the  pre-stress  on  the  specimen  through  the 
confinement.  Figure  30  shows  the  path-time  history  of  wave  and  damage  propagation.  A 
longitudinal  wave  velocity  of  10,226  m/s  has  been  determined,  which  is  in  good  agreement  with 
the  value  determined  from  the  test  with  the  spherical  projectile  (10,213  m/s).  The  best  fit  straight 
line  through  the  nucleation  sites  of  three  leading  edge  crack  centers  yields  a  damage  velocity  of 
9332  m/s.  This  is  also  very  similar  to  the  damage  velocities  obtained  with  unconfined 
specimens. 
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Figure  30.  Path-time  data  of  wave  and  fracture  propagation  in  confined  specimen, 
analysis  of  crossed  polarizers  photographs. 


The  damage  pattern  in  AlON  due  to  EOI  of  a  spherical  projectile  at  about  430  m/s  is  illustrated 
in  figure  31,  which  shows  a  selection  of  four  shadowgraphs  and  the  corresponding  crossed 
polarizers  photographs.  The  sequence  of  photographs  shows  a  fan-shaped  fracture  pattern, 
consisting  of  many  single  cracks,  emanating  from  the  center  of  impact.  The  semi-circular 
contour  indicates  that  all  cracks  are  growing  at  the  same  speed.  The  crossed  polarizers 
photographs  do  not  only  show  the  longitudinal  stress  wave  front,  but  also  the  transversal  wave 
front.  The  path-time  histories  of  the  longitudinal  and  transversal  waves  and  the  crack  front 
propagation  are  depicted  in  figure  32.  The  wave  and  crack  velocities  were  determined  by  linear 
regression  of  the  position-time  data: 

•  Longitudinal  wave  speed  (cl)  =  10,213  m/s; 

•  Transversal  wave  speed  (cj)  =  5667  m/s;  and 

•  Terminal  crack  velocity  (vcr)  =  4377  m/s. 


8.  Numerical  Simulations 


8.1  Introduction 

The  ansatz  chosen  here  for  a  numerical  modeling  of  ceramics  is  based  on  a  discretization  using 
finite  elements. 
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Figure  31.  Selection  of  four  shadowgraphs  and  corresponding  crossed  polarizers  photographs  from  baseline  tests 
with  AlON;  spherical  projectile;  test  nos.  15295-96. 


1 1 .7  |js 


Figure  32.  Path-time  history  of  wave  and  crack  propagation  in  AlON  from  baseline  tests  with 
spherical  projectile;  test  nos.  15295-96. 

In  order  to  model  the  polyhedral  micro  structure  of  the  ceramics  under  investigation  as  exactly 
as  possible,  we  first  implemented  a  code  which  is  based  on  the  theory  of  power  diagrams  (a 
generalization  of  voronoi  diagrams)  (6).  With  power  diagrams,  one  can  construct  tessellations  of 
3-D  space  which  correspond  to  a  polyhedral  cell  complex.  We  implemented  a  subsequent 
optimization  scheme  (7-9),  which  changes  the  obtained  cell  structures  according  to  the  statistics 
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of  grain  sizes  obtained  from  micrographs  in  two  dimensions,  the  only  available  source  for  the 
actual  grain  structure  of  the  materials  on  a  micrometer  scale.  The  generated  power  diagrams  are 
polyhedral  cell  complexes  in  three  dimensions  and  are  optimized  such  that  an  arbitrary  virtual 
cross  section  through  the  3-D  cell  complex  corresponds  on  average  to  the  grains’  circumference 
and  size  distribution.  This  procedure  ensures  that  our  obtained  cell,  respectively  grain  structure, 
is  as  close  to  the  experimental  evidence  as  possible. 

After  obtaining  a  3-D  grain  structure  that  corresponds  on  average  to  what  is  observed  in 
micrographs,  we  mesh  these  structures  in  three  dimensions  using  tetrahedra  in  3-D  vs.  triangles 
in  two  dimensions  (2-D).  For  2-D,  we  implemented  a  further  refinement  procedure  for  the  mesh 
which  removes  very  small  edges  and  areas  such  that  the  mesh  becomes  treatable  in  a  numerical 
FEM  ansatz.  The  influence  of  this  refinement  on  the  grains’  distribution  functions  is  negligible. 

After  having  meshed  the  cell  structures,  they  are  subject  to  FEM  simulations  using  a  commercial 
program  package  (LSDyna)  which  allows  for  a  procedure  that  can  treat  grains  as  individually 
units,  thus  allowing  an  implementation  of  failure  of  grains  (i.e.,  a  separation  of  grains). 

The  next  section  illustrates  our  ansatz  for  generating  and  optimizing  the  grain  structures  as 
polyhedral  cell  complexes  using  power  diagrams,  along  with  the  subsequent  meshing  and 
optimization  procedure.  Following  that,  the  results  of  FEM  simulations  are  displayed  and 
discussed. 

8.2  Generating  Optimized  Power  Diagrams  for  FEM  Analysis 

Figure  33  displays  our  scheme  for  obtaining  the  cell  area  and  perimeter  distribution  of  scanned 
etched  micrographs.  The  scanned  micrographs  are  manually  prepared  by  a  vector  graphics 
program  with  subsequent  analysis  using  the  KS  300  Imaging  system  by  Carl  Zeiss  Vision 
GmbH.  Figure  34  exhibits  the  only  micrographs  of  ALON  that  were  available  for  statistical 
sampling.  Here,  much  larger  samples  would  be  needed  in  the  future  to  obtain  more  reliable 
statistics. 


Vector  graphics  program  KS  300  Imaging  System, 

Carl  Zeiss  Vision  GmbH 


distribution 
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Figure  33.  Scheme  for  obtaining  the  grains’  distribution  from  etched  micrographs.  In  this 
scheme,  an  A1203  micrograph  is  shown. 


25 


2MLM 


Figure  34.  Available  etched  micrographs  of  aluminum  oxynitride  (ALON)  (7).  These  micrographs  were 
used  to  obtain  optimized  cell  complex  power  diagrams  in  three  dimensions. 


Figure  35  explains  our  optimization  scheme  for  the  3-D  power  diagram.  Several  hundreds  of 
2-D  slices  are  taken  for  each  orthogonal  direction  and  taken  into  account  for  comparison  with  the 
experimental  data.  Figure  36  exhibits  an  example  for  a  non-optimized  power  diagram,  obtained 
as  a  cut  through  the  3-D  cell  complex.  Figure  37  illustrates  in  a  float  diagram  the  algorithm  of 
the  optimization  scheme  whereas  figure  38  displays  the  efficiency  of  our  implemented 
optimization  procedure.  Figure  39  finally  exhibits  the  results  of  the  optimization  of  power 
diagrams.  Due  to  insufficient  statistical  data  on  ALON  available,  we  based  our  analysis  on 
statistical  data  of  AI2O3.  These  data  clearly  outline  the  performance  of  our  optimization  scheme. 
As  is  evident  from  the  figure,  the  optimized  structure  fits  much  better  the  underlying 
experimental  data. 


3D  model 


Figure  35.  Optimization  scheme  for  our  3-D  model  to  obtain  a  polyhedral  grain 
structure  which  corresponds  to  the  (only)  2-D  structural  information 
available  from  etched  micrographs. 
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Figure  36.  Example  of  a  generated  initial  (not  yet  optimized) 
2-D  cut  of  our  3-D  cell  complex. 


Figure  37.  Optimization  scheme  for  3-D  generated  structures.  We  use  the  first  three  central  moments  for 
comparison  of  the  generated  structures  with  the  distribution  obtained  from  etched 
micrographs. 
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Optimization  of  a  ALON  structure 


Optimization  Steps 


Figure  38.  Efficiency  of  optimization  procedure  for  ALON.  After  only  8  days  of 
optimization,  the  coefficient  ;f2 ,  which  is  a  measure  of  the  difference 
between  the  generated  structure  and  the  experimental  data,  is  for  all 
practical  purposes  very  close  to  zero. 


Figure  39.  Demonstration  of  the  result  of  the  optimization  procedure  for  the  area  and  perimeter 

distribution  for  A1203.  The  Poisson  distribution  is  the  one  obtained  from  a  non-optimized 
initial  power  diagram.  The  optimized  grain  structures  fit  much  better  the  experimental  data. 
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Figures  40  and  41  show  an  example  of  a  3-D  generated  voronoi  structure  in  the  shape  of  a  cube 
and  an  SEM  picture  of  a  typical  ceramics  surface  (in  this  case,  AI2O3,  as  no  SEM  pictures  of 
ALON  were  available).  Also  by  visual  inspection  it  is  evident  that  the  3-D  grain  structure  of  a 
ceramic  specimen  is  very  well  modelled  by  our  approach  using  power  diagrams. 


Figure  40.  Generated  and  optimized  3-D  grain  structure  of  ALON.  There 
is  a  total  of  371  gr  with  a  total  of  130,000  elements  in  three 
dimensions.  The  grains  are  initially  tied  together  by  a  contact 
force.  Dimensions  of  the  cube  are  arbitrary  units  as  the  size  of 
the  grains  can  be  scaled  to  the  desired  size  (in  this  case,  lxl 
x  1  mm). 


Figure  41.  Direct  comparison  of  a  computer-generated  granular  surface  (AlON)  and  a  SEM 
picture  of  a  typical  ceramic  grain  structure,  in  this  case  A1203. 
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Figures  42  and  43  illustrate  a  2-D  sample  along  with  the  meshing  and  the  refinement  procedure 
for  the  generated  mesh  which  leads  to  a  removal  of  short  edges  and  areas. 


Figure  42.  The  2-D  cut  through  a  3-D  cell  structure  and  subsequent  meshing  of  the  grains.  Note  that  there 
are  no  intermediate  elements  (such  as  shell  or  interface  elements)  used  in  between  the  various 
grains.  There  is  an  initial  contact  force  applied  to  the  knots  of  adjacent  interfaces  of  grains. 


Figure  43.  Illustration  of  a  mesh  optimization  procedure  in  two  dimensions  which  removes  short  edges  (and 
areas)  such  that  a  larger  time  step  in  a  subsequent  FEM  analysis  can  be  used. 
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Finally,  in  figures  44  and  45,  we  display  the  various  ways  in  which  we  are  able  to  prepare 
samples  from  our  power  diagrams’  raw  data.  Specimens  with  smooth  edges  are  possible  as  well 
as  structures  in  which  the  natural  surface  roughness  due  to  the  grains  is  kept. 


Figure  44.  Illustration  of  different  3-D  cuts  through  a  generated  power  diagram.  Structures  with 
smooth  edges  or  rough  edges  (leaving  the  surface  grains  unchanged)  are  possible. 


Figure  45.  Another  example  of  sample  preparation  from  the  raw  data  of  the  power  diagram. 

Top:  a  view  of  the  granular  structure  with  different  colors  of  grains  for  clarity. 
Bottom:  the  same  structure  but  meshed  with  tetrahedral  in  three  dimensions. 
Bottom  right:  another  enlargement  of  a  surface  area  exhibiting  the  meshed 
structure. 
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8.3  Model,  Results,  and  Discussion 


Our  approach  for  modeling  ALON  is  based  on  a  FE  discretization  and  takes  into  account  the 
microstructure  of  densely  sintered  ceramics.  The  ceramics  are  modeled  as  a  polyhedral  cell 
complex  such  that  the  statistical  properties  agree  well  with  experimental  ceramic  slices 
concerning  the  area  distribution  and  the  distribution  of  form-factors.  The  grains  of  these  cell 
complexes  were  separated  and  meshed  as  individual  objects.  Finally  these  cells  were  tied  again 
in  the  original  configuration  using  a  tiebreak  contact  implemented  in  the  commercial  FEM-code 
LS-DYNA.  Crack  opening  is  described  by  a  cohesive  softening  law  based  on  a  fracture  energy 
approach.  This  model  was  applied  to  an  8-  x  8 -mm  AlON-ceramics  2-D  plain  strain  model  hit 
by  a  steel-object  at  high  velocity  (430  m/s).  A  simple  isotropic,  linear  elastic  material  model  was 
used  for  the  grains  and  for  Young’s  modulus  E.  The  maximum  stress  upon  which  intra-granular 
failure  occurs  and  the  fracture  energy  of  the  contacts  were  taken  from  macroscopic  data.  The 
results  of  the  simulation  are  shown  in  figure  46. 


Figure  46.  Left:  AlON  plate  LS-DYNA  model,  8x8  mm,  0.3  ps  after  being  hit  by  a  steel 
sphere  at  430  m/s.  Right:  same  model  after  0.6  ps;  the  colors  mark  the  pressure 
level  -1  GPa  to  +1  Gpa;  the  blue  line  marks  a  zone  where  most  contacts  have 
already  failed. 

The  zone  where  fully  developed  cracks  between  individual  grains  can  be  observed  is  marked  by 
a  blue  line.  The  growth  rate  of  this  zone  is  about  half  the  velocity  of  the  elastic  wave  speed, 
which  compares  well  with  experimental  observations  (figure  31).  Under  the  assumption  of  a 
constant  wave-velocity  and  symmetric  behavior,  one  can  qualitatively  compare  the  model  at 
0.6  ps  with  the  experiment  at  6.7  ps. 

Note  that  our  model  does  not  include  any  arbitrary  additional  elements  at  the  grain  boundaries 
such  as  interface  elements  or  shell  elements,  etc.  In  this  case,  one  would  always  change  the  ratio 
of  volume  to  surface  within  the  grains  and  always  overestimate  the  influence  of  the  grain 
boundary.  The  grain  boundaries  are  just  areas  of  local  atomic  disorder  of  adjoining  single 
crystals.  The  size  of  this  area  is  negligibly  small  (a  few  Angstroms)  compared  to  the  typical  size 
of  a  grain.  Therefore  we  use  an  ansatz  with  contact  forces  at  grain  surfaces  which  initially  glue 
the  grains’  knots  at  surfaces  together  and  which  does  not  use  any  additional  and  artificial 
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interface  elements.  As  we  have  no  experimental  data  for  this  crystal  contact  force,  we  used  and 
tested  different  material  models  implemented  in  the  commercial  code  LS-DYNA. 

In  a  first  attempt,  we  used  material  no.  96  in  the  commercial  code  LS-DYNA  which  allegedly  is 
a  model  suited  for  the  description  of  brittle  damage.  This  model  requires  a  number  of  material 
parameters  as  input,  namely: 

•  Density 

•  Young’s  modulus  E 

•  Poisson  number  n 

•  tensile  and  shear  stress  limit 

•  fracture  toughness 

•  compressive  yield  stress 

We  only  had  available  E  and  the  fracture  toughness  for  AlON.  The  tensile  and  shear  limit  were 
taken  as  8T05  Pa  and  2-105  Pa,  respectively,  which  is  a  ratio  of  four  in  reference  to  concrete. 

The  compressive  yield  stress  was  estimated  as  being  close  to  that  of  AI2O3.  We  tested  this  model 
in  a  total  of  19  different  parameter  variations  but  in  all  variations  it  led  to  strong  mesh  distortion 
and  subsequent  termination  of  simulation  due  to  numerical  instabilities.  Therefore,  we  finally 
used  a  simple  piecewise  elastic  plastic  material  model  for  AlON  provided  by  the  LSDyna  library 
using  the  material  parameters  of  ALON  (p  jli,  E,  and  at  [  /  0] )  for  the  impacted  specimen,  and 
steel  (from  ww.matweb.com)  for  the  impactor.  We  then  observe  the  correct  wave  velocity  in  the 
material  (-10,040  m/s)  and  a  velocity  of  propagating  cracks  of  -4530  m/s  which  is  in  good 
agreement  with  the  experiments  (see  figures  3 1  and  32). 

Note  that  only  by  using  this  specific  way  of  modeling,  namely  using  tiebreak  contacts  between 
the  grains,  we  are  able  to  really  model  failure  of  the  material  and  to  avoid  any  artificial  numerical 
behavior,  as  can  be  seen  when  modeling  the  target  as  a  continuum  without  the  microstructural 
details. 

As  an  illustration  of  this,  we  simulated  a  full  continuum  model  that  does  not  contain  any 
microstructural  details  using  the  commercial  package  Autodyn,  and  show  that  this  model 
depends  strongly  on  the  chosen  mesh  resolution  (see  figure  47).  In  the  left  picture,  the  mesh 
resolution  was  1  mm;  in  the  right  one  it  was  0.5  mm.  The  results  of  failure  are  completely 
different  in  these  two  pictures  with  otherwise  identical  simulation  conditions.  Also,  the 
simulation  using  smooth  particles  does  not  lead  to  any  improvement  in  this  respect  (see  figure  48). 
Here,  there  are  even  more  artifacts  visible  in  respect  to  failed  elements. 
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Figure  47.  Insufficiency  of  generic  FEM  approaches.  A  steel  impactor  hits  an  AlON  tile. 
Red:  failed  elements.  Starting  configuration  of  an  EOI  at  vP  =  50  m/s. 


Figure  48.  AlON  modeled  with  SPFI  particles  and  a  resolution 
of  1  mm.  The  result  is  different  from  the  one  in 
figure  47  (left),  although  all  simulation  conditions 
are  the  same. 


Figures  49-5 1  show  a  series  of  snapshots  from  an  impact  simulation  with  a  2-D  system  of  size 
9x9  mm  at  vp  =  430  m/s.  One  can  see  nicely  the  onset  of  developing  cracks  in  the  material 
until  complete  failure  with  several  chunks  of  failed  material.  Finally,  in  figure  52,  the  actual  size 
limitations  in  the  simulations  are  shown.  Only  a  very  small  area  of  the  ceramic  tile  can  be 
simulated  when  the  microstructure  is  taken  into  account  explicitly. 
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Figure  49.  Impacted  AlON  at  450  m/s  with  a  steel 
impactor,  1.8  10“7s  after  impact. 


Figure  50.  Impacted  AlON  at  450  m/s  with  a  steel 
impactor,  7.4  10“7  s  after  impact. 
Comminution  of  the  impact  zone  has 
already  started. 


8.4  Outlook  (Numerical  Simulation) 

A  fundamental  problem  which  we  are  facing  with  our  FEM  approach  is  the  fact  that  the 
algorithm  which  is  implemented  for  the  contact  force  is  extremely  expensive,  thus  rendering 
simulations  of  3-D  systems  or  intra- grain  failure  practically  impossible.  Even  for  2-D  systems, 
we  are  facing  extremely  small  time  steps  and  long  simulation  runs  with  more  than  500,000 
elements  along  with  many  numerical  stability  problems  due  to  strong  mesh  distortions. 
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Figure  5 1 .  Impacted  A10N  at  450  m/s  with  a  steel 
impactor,  4.9  10”6s  after  impact. 
Failure  of  the  ceramic  with  different 
fracture  pattern  propagating  through 
the  ceramic  tile. 


Figure  52.  Snapshot  of  an  impact  experiment  with  a  spherical  steel  impactor  and  the  magnified  simulated  2-D 
sections  of  the  AlON  ceramic  tile  which  exhibit  the  damage  zone  after  the  impact.  The  simulation 
model  is  a  7.8-  x  7. 8 -mm2  section  of  the  original  100-  x  100-mm2  tile.  Using  an  elastic  plastic  material 
model,  the  wave  velocities  in  the  experiment  and  in  the  simulations  are  the  same.  The  simulation 
snapshot  is  taken  after  the  shock  wave  has  crossed  about  70%  of  the  simulated  part  of  the  tile. 
Displayed  minimum  and  maximum  pressure  is  ±1  GPa. 


What  is  necessary  in  these  endeavors  is  the  development  of  an  apt  and  more  efficient  code 
structure,  optimized  for  simulating  failure  in  a  mesh  without  using  additional  interface  elements. 
Also,  more  etched  2-D  viewgraphs  displaying  the  granular  structure  of  AlON  would  be  very 
helpful  for  an  enhanced  optimization  of  the  power  diagrams.  A  different  method  for  simulating 
facture  and  failure  in  ceramics,  developed  by  Steinhauser  et  al.  ( 11-13 ),  based  on  overlapping 
particles  interacting  via  potentials,  shows  very  promising  results  for  further  investigations  along 
this  direction.  However,  with  the  limitations  of  time  and  resources  within  this  project,  further 
development  of  this  method  was  not  feasible. 
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9.  Conclusion 


The  EOI  test  method  was  applied  in  order  to  visualize  wave  and  damage  propagation  in  materials 
for  transparent  armor.  The  influence  of  bonding  layer  thickness  on  damage  evolution  in 
Starphire  glass  laminates  was  examined.  The  high  resolution  of  the  high-speed  photographs 
allowed  for  the  determination  of  the  stress  wave  time  delay  during  the  transition  through  the 
bonding  layers.  It  is  expected  that  the  capabilities  of  the  experimental  method  help  with  the 
development  of  damage  models  and  that  the  combination  of  experimental  and  computational 
modelling  results  can  eventually  guide  materials  and  laminates  design. 
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Appendix.  Complete  Sets  of  High-Speed  Photographs  From 
Edge-on-Impact  (EOI)  Tests 
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Figure  A-l.  Complete  set  of  high-speed  photographs  from  EOI  test  of  laminated  Starphire  specimen  with 
saw  tooth  interface,  shadowgraph  configuration;  test  no.  15724. 
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Figure  A-2.  Complete  set  of  high-speed  photographs  from  EOI  test  of  laminated  Starphire  specimen  with 
saw  tooth  interface,  crossed  polarizers  configuration;  test  no.  15725. 
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Figure  A-3.  Complete  set  of  high-speed  photographs  from  EOI  test  of  laminated  Starphire  specimen  with 
corrugated  interface,  shadowgraph  configuration;  test  no.  15727. 
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Figure  A-4.  Complete  set  of  high-speed  photographs  from  EOI  test  of  laminated  Starphire  specimen  with 
corrugated  interface,  crossed  polarizers  configuration;  test  no.  15728. 
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Figure  A-5.  Complete  set  of  high-speed  photographs  from  EOI  test  of  laminated  Starphire  specimen  with  wave 
shaped  interface,  shadowgraph  configuration;  test  no.  15731. 
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Figure  A-6.  Complete  set  of  high-speed  photographs  from  EOI  test  of  laminated  Starphire  specimen  with 
wave  shaped  interface,  crossed  polarizers  configuration;  test  no.  15730. 
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